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Separate spin evolution quantum hydrodynamics is generalized to include the Coulomb exchange 
interaction. The Coulomb exchange interaction is considered as the interaction between the spin- 
down electrons being in the quantum states occupied by one electron, giving main contribution in the 
equilibrium. The generalized model is applied to study the non-linear spin-electron acoustic waves. 
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I. INTRODUCTION 


Spin evolution in quantum plasmas has been consid¬ 
ered for a long time. Fundamental equations of many- 
particle smn.-^2 quantum hydrodynamics were derived 
in 2001 [l|, [3- They found many applications at the 
study of waves and instabilities in spin-1/2 quantum plas¬ 
mas. Among them we find the appearance of the spin- 
plasma waves with frequencies near the cyclotron fre¬ 
quency Q, 0, [il, [6[. Later it was demonstrated that 
the quantum Bohm potential existing in the magnetic 
moment evolution ^nation shifts the frequency of the 
spin-plasma waves Q, [11 • This shift is proportional to 
the square of the wave vector. The annihilation interac¬ 
tion between electrons and positrons gives a shift of the 
spin-plasma wave frequency on a constant proportional 
to the magnetic moment density of the spin-1/2 quan¬ 
tum electron-positron plasmas [H. Generation of waves 
in magnetized plasmas by neutron beams via the spin- 
spin and spin-current interactions of the neutron spins 
with the spins and electric currents of the electrons and 
ions of the plasmas was considered in Ref. [bl . Spin parts 
of the hydrodynamic vorticity and helicity of the spin-1/2 
quantum plasmas were derived in Refs, [ol, [Io|. Conser¬ 
vation of the full helicity at the charge-charge and spin- 
spin interactions was demonstrated there. Spin parts of 
the vorticity and helicity for spin-1/2 electron-positron 
plasmas were obtained in Ref. [H . These phenomena ap¬ 
pear along with the change of properties of well-known 
plasma phenomena (see for instance Refs. [HI, 0, El, 

El). 

Different methods of derivation [HI, [ibl and gener¬ 
alization 0, 0, 0, 0 of the spin-1/2 quan¬ 
tum hydrodynamics were presented in literature. Kinetic 
models for the spin-1/2 quant urn plasmas were presented 
as well 0, [HI (see also reviews [Ij, [22|). All these mod- 
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els consider electrons as a single fluid. It corresponds to 
the multi-fluid plasmas, where each species is considered 
as a fluid. The Pauli equation allows to find different 
form for the spin-1/2 quantum hydrodynamics, where 
we have two fluids of electrons: the spin-np electrons and 
the spin-down electrons [23| , 0 ^ EH • Early papers [23| , 
0 did not show full picture of electron evolution and 
partially mistreat coefficients in the spin-spin interaction 
force. The separate spin evolution of the electrons, in ac¬ 
cordance with the Pauli equation, was obtained in Ref. 


It was demonstrated in Ref. [H| that the Fermi pres¬ 
sures for the spin-up electrons and the spin-down elec¬ 
trons are different. It leads to existence of new phenom¬ 
ena. The spin-electron acoustic wave was found in Ref. 
[HI at wave propagation parallel to the external magnetic 
field. Oblique propagation of the longitudinal waves at 
the separate evolution of the spin-up electrons and the 
spin-down electrons was considered in Ref. [26|- Ex¬ 
istence of two kinds of the spin-electron acoustic waves 
(SEAWs) was demonstrated in this regime. Properties of 
the SEAWs in two dimensional structures were studied in 
Ref. [23 • The SEAWs in the two-dimensional plane-like 
electron gas and the electron ga s on the cylindric sur¬ 
face were considered in Ref. [^. Kinetic model of the 
SEAWs was considered in Ref. [Hl^ where the Landau 
damping of the SEAWs was calculated. It was demon¬ 
strated that the Landau damping of the SEAWs is small. 
Therefore, the SEAWs are slowly damping waves. 


Spins of electrons affect the plasma dynamics even if 
we do not consider spin evolution. It is enough to include 
distribution of electrons on different spin states to find a 
change in the equation of state. The distribution of elec¬ 
trons on spin states affects the Coulomb exchange inter¬ 
action as well. The exchange interaction was considered 
at the first steps of the development of the many-particle 
spin-1/2 quantum hydrodynamics The exchange in¬ 

teraction attracts a lot of attention in recent research (see 
for instance Refs. 0-0). This research are grounded 
on the long experience of the Coulomb exchange interac- 
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tion study [35|-[4l|, along with recent applications of the 
exchange interaction to different plasma phenomena (see 

for instance M, Ii3)- 

In effort to study the influence of the Coulomb ex¬ 
change interaction on the properties of the spin-electron 
acoustic waves we develop generalization of the separate 
spin evolution quantum hydrodynamics (SSE-QHDs) [25| 
containing contribution of the exchange interaction. As 
in Ref. we focus our attention on degenerate electron 
gas. 

In this paper we consider the Coulomb exchange inter¬ 
action, which arises as the interaction between the spin- 
down electrons being in the quantum states occupied by 
one electron. It is considered in Ref. in terms of 
single fluid model of electrons. One of the features of 
this paper is the application of the exchange interaction 
obtained in Ref. Q to two-fluid model of electrons [25| . 

We apply the developed model to the non-linear 
SEAWs, particularly to the soliton formation. 

Non-linear waves related to the ion-acoustic waves are 
still under consideration [ij, while this paper is dedi¬ 
cated to non-linear waves related to the recently found 
spin-electron acoustic waves [25| . 

This paper is organized as follows. In Sec. II we 
present the QHD model with separated spin-up elec¬ 
trons and spin-down electrons containing the Coulomb 
exchange interaction. In Sec. Ill we describe method of 
derivation of the spin-electron acoustic soliton from the 
developed in Sec. II model. In Sec. IV we present anal¬ 
ysis of properties of the spin-electron acoustic soliton in 
quantum plasmas with no account of the exchange inter¬ 
action. In Sec. V we describe the spin-electron acoustic 
soliton with the account of the exchange interaction. In 
Sec. VI brief summary of obtained results is presented. 


II. MODEL 


The Pauli equation is a set of two equations describ¬ 
ing evolution of two wave functions, one is for spin-up 
state of electron and another one is for spin-down state of 
electron. Therefore, the evolution of system of electrons 
can be described in terms of two-fluid model of electrons 
with different spin projection. This is called the separate 
spin evolution quantum hydrodynamics 1^. . Corre¬ 

sponding kinetic model is obtained as well . However 
these models are derived in the self-consistent field ap¬ 
proximation. In this paper we make the next step in the 
development of the SSE-QHD. We include the Coulomb 
exchange interaction in the SSE-QHD. 

The time evolution of the concentrations of the spin-up 
electrons and the spin-down electrons obeys the continu¬ 
ity equations with nonzero right-hand side 


dtn^ + V(ntV|) = ^{S^By - SyB^), 

and 

dtn^ + VKv^) = -^{S,By - SyB,), 


( 1 ) 


where and ni (v^ and v^) are the particle concen¬ 
trations (the velocity fields) of the spin-up electrons and 
the spin-down electrons, fie = —9^^^ is the magnetic mo¬ 
ment of electron, and g = 1 a/liir) = 1.00116, where 
a = 1/137 is the fine structure constant, which gets into 
account the anomaious magnetic moment of eiectron, e 
(rUe) is the eiectron charge (mass), Ti is the Planck con¬ 
stant, c is the speed of light, B = {Bx^By^Bz} is the 
magnetic field. The particle concentrations appear as 
the quantum mechanical average of the corresponding 
wave functions, which are the elements of the Pauli spinor 
wave function, Ug = ('0*'0s), with s =t or j.. These con¬ 
centrations are related to the spin-up electrons and the 
spin-down electrons separately. It appears in the accor¬ 
dance with the spinor structure of the Pauli equation, 
which governs the evolution of the spinor wave function 

7 / = ( ^^ ). Considering the evolution of the upper and 

lower elements separately we find the separate descrip¬ 
tion of the spin-up eiectrons and the spin-down eiectrons. 
The veiocity fieids and appear at the averaging 
of the corresponding operators with 7 /^ and 7 /^: = 

( 1 /rig)(( 7 /*P '05 -h c.c.)/2me), where c.c. stands for the 
compiex conjugation. Equations m and ([ 2 ]) contain pro¬ 
jections of the spin density Sx and Sy. Each projection of 
the spin density is defined as a mixture of the spin-up and 
spin-down wave functions: Sx = '0*cra^'0 = 01 
and Sy = ^(010-1^ — 0|0^). Therefore, these 

quantities are not reiated to different species of eiectrons 
having different spin direction. Sx and Sy describe simui- 
taneous evoiution of both species of eiectrons. Equations 
of evoiution of Sx and Sy were derived in Ref. as a 
part of the set of SSE-QHD equations. We do not study 
the spin evoiution in this paper, so we do not describe 
equations for Sx and Sy^ which can be found in Refs. [25[ 
and [ 2 ^. 

Erom the continuity equations o and we see that 
the numbers of eiectrons in each subspecies can change 
due to the spin-spin interaction and the interaction of 
spins with the externai magnetic fieid. However, the fuii 
number of eiectrons Ue = conserves in this modei. 

In this modei we have two Euier equations. We use 
the subindex s or 0 to present them as one equation 


, .7 (Vn,)2 

mns[dt + + Vps - I — - 


= geUgfE-l- -[v5,B] ) -^Fex, 


-h ^{SxS/ Bx + SyS/By) 




ififle 


^ i^(M)xBy J (^M)yBx) ^ ^ {SxBy SyBx), 

(3) 


( 2 ) 
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where in coefficients ± and ^ we have the upper sign 
for the spin-up electrons and the lower one for the spin- 
down electrons. In formula dU we use qe = —e for the 
electron charge, Ps for the pressure of the spin-up and 
spin-down electrons. We also apply J(m)x ^iM)y 
for the elements of the spin current tensor . Vectors 
J(M)cc arid ^{M)y have the following explicit forms 


^(M)x - ^(Vt + ^ 




(4) 


and 




The Euler equations ([3]) describe the momentum evo¬ 
lution complicated by the unconservation of the numbers 
of the spin-up and spin-down electrons. 

We describe now the physical meaning of different 
terms in the Euler equations ©• The first term is the 
continual derivative of the velocity field v^. The second 
term is the gradient of pressure. The explicit form of the 
pressure we present and discuss below. The third term 
is proportional to the square of the Planck constant. It 
is a combination of the spatial derivatives of the particle 
concentration Ug up to the third derivative Arig. This 
term is called the quantum Bohm potential. It is related 
to the wave nature of the electron. 

On the right-hand side of the Euler equations m we 
present the force fields of different nature. The first term 
presents the Lorentz force describing the interaction of 
charges with the electromagnetic fields. The electric field 
contains two parts E = —V0 — 9tA/c, where (j) and A 
are the scalar and vector potentials of the electromag¬ 
netic field. The first of them is the potential part giving 
contribution in the longitudinal waves, which we consider 
in this paper, and the second one is the vortical part. The 
second part of the Lorentz force describes the interaction 
of the moving charges with the magnetic field. The sec¬ 
ond term on the right-hand side of the Euler equations 
([3|) is the exchange part of the Coulomb interaction. The 
contribution of the Coulomb exchange interaction in the 
separate spin evolution QHD is in the center of attention 
of this paper. We discuss its explicit form below. 

The third and fourth terms are the parts of the spin- 
spin interaction force. The magnetic moments related 
to the spin of electrons create the magnetic field. This 
magnetic field acts on the magnetic moments of other 
electrons and leads to existence of this force in the Eu¬ 
ler equation. In terms of the SSE-QHD this force field 
splits on two terms. The first (the second) of them de¬ 
scribes the interaction of the z projection (the x- and y 
projections) of the magnetic moments with the nonuni¬ 
form z projection (the x- and y projections) of the mag¬ 
netic field. The last two terms on the right-hand side of 
the Euler equations (|3|) are related to the unconservation 
of the numbers of the spin-up and spin-down electrons. 
The first of them appears at the derivation of the Euler 


equation for the momentum density The second of 

them arises at the application of the continuity equation 
during the extraction of dtVg from dt{nsVs)- Hence, it 
is proportional to the right-hand side of corresponding 
continuity equation. 

Study of the Coulomb exchange interaction in the elec¬ 
tron gas has a long history. Recently, it was shown that 
the exchange interaction force field stron gly depends on 
the spin polarization of the electron gas [^. Distribu¬ 
tion of the partially spin polarized electrons is depicted 
in Eig. [H This distribution splits the electrons on differ¬ 
ent groups, which demonstrate different exchange inter¬ 
actions. Eig. ([T]a) shows two pairs of electrons. We see 
a spin-down electron being in a state with energy E G 
{^Fe{up)^^Fe{down)] interacting with two electrons hav¬ 
ing different spin direction and being in the same quan¬ 
tum state with energy E' G [^^£Fe{up)\i where £Fe{up) = 

( 67 r 2 not)^/^?i^/ 2 m, eFe(down) = I2m. The 

strengths of these interactions are the same, but they 
have opposite signs. Hence, it gives zero contribution in 
the force field. 

In Eig. (Pb) we have a similar situation, but a chosen 
spin-up electron is in a quantum state with energy E G 
[0,£Fe{up)]' H gives zeroth contribution in the force field 
either. If we consider a spin-down electron in a quantum 
state with energy E G [0, £Fe{up)] we find the same result. 

Fig. (Pc) shows the regime giving non zero contribu¬ 
tion in the force field. The regime presented in Eig. P 
c) was considered in Ref. [s^. Here, we have interaction 
of two spin-down electrons being in quantum states with 
energies E G (^Fe(np),^Fe(do^/;n)]- In this regime the spa¬ 
tial part of the wave function is antisymmetric relatively 
permutation of these two particles. We have same sign 
of interaction for all such pairs of electrons. It is impor¬ 
tant to underline that this regime involves the interaction 
of spin-down electrons only. Hence we should substitute 
this force field, found in Ref. [s^, in the Euler equation 
for the spin-down electrons. This regime gives no contri¬ 
bution in the Euler equation for the spin-up electrons. 

The regime of the electron interaction depicted in Eig. 
P d) describes the electrons having opposite spins and 
located in the same quantum state. In this case the 
spatial part of the wave function is symmetric relatively 
to the permutation of two particles even without addi¬ 
tional symmetrization V^p.p. (ri, r 2 ) = V^p. (ri)' 0 p. (r 2 ) = 
T/^p.p. (ri, r2). Therefore, there is no exchange interaction 
in this regime. 

The force of the exchange interaction of the spin-down 
electrons being in quantum states occupied by one elec¬ 
tron was found in Ref. [s^. The result was presented 
in terms of the concentration of all electrons, while it 
involves the interaction of spin-down electrons only. To 
substitute this force in the Euler equation for the spin- 
down electrons it is necessary to rewrite the force field 
in terms of the spin-down electron concentration. In the 
equilibrium state we have the following relations between 
the spin polarization 77 , the spin-up electron concentra¬ 
tion the spin-down electron concentration and 
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FIG. 1: (Color online) The figure shows the different combi¬ 
nations of pair of the spin-up and spin-down electrons in the 
degenerate electron gas giving different contributions in the 
collective effect of the Coulomb exchange interaction. In this 
figure we apply the Fermi energies of the spin-up and spin- 
down electrons: epeiup) = /2m, Speidown) = 


the full concentration of electrons Ue'. = (1 — rj)ne/2, 
ni = (1 + 7^)ne/2. Applying these formulae we can make 
the required representation of the exchange interaction 
force: 


'Pex,U = ( 6 ) 


where 


C3D = (i + r/F"-(i-r/F", (7) 


and 


X = CsD 


23 


TT (1 + 77 )^/^ 


= 23 


3 / 3 / (l-r?)V3 \ 

V TT V (1 + 77 )^/^/ 


( 8 ) 


The electromagnetic field presented in the hydrody¬ 
namic equations satisfy the Maxwell equations: 


waves k || we find from equation ([TT]) that the per¬ 
turbation of the magnetic field is equal to zero = 0. 
Hence the set of SSE-QHD equations ([I])-(|5j) simplifies to 
the following set of equations. 

We consider interval of the equilibrium concentrations 
from 77-0 = 10^^ cm“^ to no = 10^^ cm“^. We drop 
the contribution of the quantum Bohm potential, which 
reveals itself at larger concentrations. 

In this regime we have conservation of the electron 
number for both subspecies 



dtn^ + V(n-|-V|) = 0, 

(13) 

and 

dfUi + V(n4,V4,) = 0. 

(14) 

The conservation follows from the zeroth value on the 
right-hand sides of the continuity equations ([T3)) and 

da- 

Simplified Euler equations appear as follows 

mn^{dt + v^V)v| + 



= ge%E+ —ni-[vi.,B], 

(15) 


and 


mni{dt + v;V)v; + 

= geniE + X9ev/^Vn4, + ^n4,[v4,,B], (16) 

where we present the explicit form of the Coulomb ex¬ 
change interaction. The magnetic field B in equations 
m and ([T6|) is the external magnetic field. 

The electric field in simplified Euler equations ([T5j) and 
([T6|) is the quasi-static electric field. Hence it obeys the 
Poisson equation 


VE = 47r(erii - eriet - euei), (9) 


VE = 47r(enio - eue^ - euei), (17) 


VB = 0, 

(10) 

VxE = -IatB, 
c 

(11) 


and 


and the eddy-free condition 

VxE = 0. (18) 

In equilibrium state we have Uio = noe = 77-ot + • 

Considering quantum spin-1/2 plasmas researchers 
usually apply equation of state for unpolarized electrons 


V X B = -dtE 
c 


H-Tl (^“^“tVat + 9anaiVo4,)+47r ^ VxMa, (12) 

a=e,i a=e,i 

where = {lJ.aSax, l^aSay, l^a{nat-nai)} is the magne- 
tization of electrons in terms of hydrodynamic variables. 

We consider the wave propagation parallel to the ex¬ 
ternal magnetic field. Since we consider the longitudinal 


(37r^)3 5 

Punpol ~ Z 77-3^ 

5 m 


(19) 


see Ref. ji^ . 

One can include the contribution of the spin polariza¬ 
tion of the degenerate electron gas in the pressure in the 
single fluid model of the three dimensional electron gas 

Q, 


Pe = 1 ^ 30 (StT^) 3 


0 ^ 
h^ni 

brue 


( 20 ) 
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where 


= + + ( 21 ) 

is the coefficient describing the spin polarization of the 
electron gas. 

However, we deal with the separate evolution of the 
spin-up electrons and the spin-down electrons. There¬ 
fore, we apply the partial pressure caused by each species 
of electrons Hence, we have 


Pet = (67r2)" 



brrie ’ 


for the spin-up electrons, and 


Pel = (OtT^)" 



brrie ’ 


( 22 ) 


(23) 


for the spin-down electrons. Sum of Pet Pei gives us 
Pe presented by formula (|2Q]) if we include the relation 
between concentrations = (1 — r])nel2 and = (1 -h 
■n)nel2. 

Difference of the spin-up electron concentration and 
the spin-down electron concentration An = not ~ 
is caused by the external magnetic field. Since elec¬ 
trons are negatively charged their spins have prefer¬ 
able direction opposite to the external magnetic field 
= tanh(i^) = -tanh(%^), t] =\ An \ 

where j2m is the Fermi tempera¬ 

ture in units of energy, so we do not write the Boltzmann 
constant. 


III. PERTURBATION EVOLUTION 


Following Ref. [13 we introduce an expansion of the 
hydrodynamic parameters on small parameter 5 

Us = nos + enxs + (26) 

— 0 “ 1 “ ^ '^2sz') ( 27 ) 

and 

0 = 0 + 501 + 5^02, (28) 

where 0 is the potential of the electric field E = — V0. 

We substitute formulae ([26|) - ([28|) in equations m- 
(pT|) In the leading order on the small parameter 5 we 
find the following relations between the perturbations of 
the particle concentrations, the velocity fields, and the 
potential of the electric field 


V 

— nis, 

nos 

(29) 

—enot0i 

(30) 

—eno^0i 

“ me(y2 - C/2)’ 

(31) 


where we have applied the following notations 


and 


We consider the propagation of the non-linear pertur¬ 
bations parallel to the external magnetic field B = BqBz. 
Here, the plane wave soliton propagates parallel to the ex¬ 
ternal magnetic field. Hence it parameters depend on the 
single coordinate and we have one dimensional perturba¬ 
tion. We focus our attention on the non-linear waves 
related to the SEAWs. We consider ions as the motion¬ 
less positively charged background. To find the soliton 
solution we apply the perturbation technic developed by 
Washimi and Taniuti in Ref. [13 • This technic is widely 
applied in recent research of wave phenomena in plasma 
physics. In this paper we apply it to find the spin-electron 
acoustic soliton. 

In equations ([T3|) - ([T6|) we make transition to variables 
^ and r defined as follows 

^ = ei{z-Vt), (24) 

and 

T = (25) 

where 5 ^ 1 is a dimensionless parameter. 


Uf = 


3m2 


(67r^no^)3 - —n, 

me 


0^- 


(33) 


Below, we find that V‘^ — U^ and V‘^ — U^ have different 
signs. Hence perturbations ni^ and ni^ have different 
signs either. 

The Poisson equation in the leading order on 5 appears 
as 


= 0- (34) 

Substituting formulae and m in equation (|34|) 
we find the velocity of the perturbation propagation in¬ 
troduced in formula (|24)) : 

y" = —{no^uf + noiU^). (35) 

noe 

We need to find the explicit form of perturbations ni^, 
nil, ^^4 01. To this end, we consider the hydrodynamic 
equations ([T3]) - ([T7|) in the next order on small parameter 
5. In this regime we find n 2 s, V 2 sz in terms of 0i. We 
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FIG. 2: (Color online) The figure presents the spin polariza¬ 
tion dependence of the soliton profile of the spin-down elec¬ 
tron concentration. We present the soliton prohle in terms 
of 'Ed = {Tiii/noe)/{V^Uo/vFe), whcrc nil appears at the 
substitution of solution (l38l) in formula (|3T]). 
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FIG. 3: (Color online) The figure presents the spin polariza¬ 
tion dependence of the soliton profile of the spin-up electron 
concentration. Similarly to the previous figure, we present 
the soliton profile in terms of Eu = (nii/noe)/(V^Uo/ vfc), 
where nif appears at the substitution of solution (l38l) in for¬ 
mula (l3Q|) . 



FIG. 4: (Color online) The figure shows the spin dependence 
of the square of the dimensionless perturbation velocity in 
self-consistent field approximation defined by formula (1391) . 


substitute expressions for the second order perturbations 
of the particle concentration n 2 s in the Poisson equation 

= 47re(n2t + ^ 24 ,) (36) 

obtained from dm in the second order on 5 . 

Finally we find equation for the first order perturbation 
of the electric field potential 

Adr4>i + dlh - Bd^ll = 0- (37) 

It is the KdV equation. Coefficients A and B are rather 
huge. We present and discuss them in Secs. IV and V. 
Let us mention that coefficient A is positive ^4 > 0. 

Applying the substitution ^ = ^ — Uqt we find a soliton 
solution of KdV equation (|37|) 

3AUo 1 

2B cosh2(iy]4C7oO' 


Parameter Vo is the velocity of the soliton propagation. 
From formula (|38|) we see that the amplitude of the soli¬ 
ton is proportional to the soliton velocity Vo, while the 
soliton width A is proportional to the inverse square root 
of the soliton velocity Vo- 

We do not consider the ion contribution assuming ions 
as motionless. Therefore we have a limit on the perturba¬ 
tion velocity V.It should be larger than the sound veloc- 


ity Vs'. V :$>Vs = where vpe = {Sir^noy^^h/m 

is the Fermi velocity. Fig. 0] shows that at large spin 
polarization r] ^ 1 the perturbation velocity becomes 
rather small vpe ^ V. In this regime the ion motion is 
essential. Therefore, we do not consider r] > 0.9 at the 
soliton study. 

To find the spin electron acoustic soliton we apply a 
small perturbation evolution method. Consequently we 
have condition for the soliton amplitudes nis rios- this 
condition gives restrictions for the velocity of the soliton 
propagation. For instance, considering the spin-up elec¬ 
trons we obtain ^ ^ Vo ^ 




Aeji 


IV. SPIN ELECTRON ACOUSTIC SOLITON IN 

THE SELF-CONSISTENT FIELD 
APPROXIMATION 

In this section we present the analysis of the existence 
and properties of the spin-electron acoustic soliton. To 
make this analysis simpler we consider the self-consistent 
field approximation. That means we drop the exchange 
interaction till the next section. 

If we do not consider the exchange interaction we 
should drop the second term in the definition of V^ (see 
formula (|33|) ). 

It is useful to represent the velocity of perturbation 

V, given by formula (|35]) . in terms of the Fermi velocity. 
Hence, in the self-consistent field approximation, we find 

• w‘^ where 

= ^{l-ri)i{l + ri)i[{l - 7])^ + (l + r?)3], (39) 
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Coefficient w‘^ describes the dependence of the perturba¬ 
tion velocity on the spin state of the electron gas. 

Coefficients A and B in the KdV equation (|37|) , in the 
self-consistent field approximation, appear as follows 

I 1 + ^ 1 

I [w^ “ (1 “ “ (1 + ) 


Vs f _1_ 

^Fer^e \ (1 - + (1 + - 1 ]^ 


+ 


and 


(l + r?)3[i(l-r?)l[(l-7?)3 +(l + r;)3] - 1]2 


(40) 


B = 


File 


f hV-v) 

I 777 ^ + 1 ( 1 - 77)3 


2 77 ;^ — (I — 77 )^ 


[w‘^ — (1 + 77 ) 


1 

2 


+ 5 ( 1 + 77)3 1 

— (1 -h 77)^ J ^ 


(41) 


where roe = vpe!^ is the Debay radius. 

Coefficients A and B can be written as ^4 = Aq • ^( 77 ) 
and B = Bq • B{r]), where Aq = and Bq = 

3«Fe^Dee/TOe- 

We see that coefficient A does not depend on the equi¬ 
librium particle concentration. 

The spin-electron acoustic soliton exists at the inter¬ 
mediate spin polarizations. It disappears at 77 ^ 0 and 
77 ^ 1. Formally it corresponds to ^4 ^ oc at 77 ^ 0 and 
77 ^ 1 . Coefficient B becomes infinite 5 ^ 00 at 77 0 

and 77 ^ 1 either. 

At the numerical analysis we consider the spin-electron 
acoustic soliton as the perturbations of the particle con¬ 
centrations of the spin-up electrons and the spin-down 
electrons substituting solution (l38]) in formulae (|3Q|) and 

(|3T]). 

The spin-electron acoustic soliton, for in¬ 
stance, for the spin-down electrons, arises 
as Tiid = nodV-5A/B){VSUo/vFe)[w^ - (1 + 

77 )^/ 3]“1 cosh-^(c 7 / 2 A), where A = VFe/VSUoArDe is 

the soliton width. 

Relative perturbations of the particle concentrations 
are proportional to the velocity of the soliton 
propagation Uq in units of the Fermi velocity vpe^ Pre¬ 
senting soliton profiles in figures we apply the effective 
amplitude = {nis/noe)/{VSUo/vFe)- 

Fig. dD) (Fig. 0) shows that the particle concentra- 
tion perturbations are negative (positive) for the spin- 
down (spin-up) electrons. Applying notations accepted 
in some areas of the condensed matter physics and optics 


we can call the negative (positive) solitonic perturbations 
as the dark (bright) soliton. 

For the particle concentration noe = 10^^ cm“^ we 
have VFetV^ = 7.9 x 10^ cm/s. Choosing Uq = 10“^7;Fe 
we obtain 77 is = 0.5(1 ± 77 ) 10 ^^Ss, where Es is shown in 
Fig. (j2j) for the spin-down electrons and in Fig. (j3]) for 
the spin-up electrons. Factor 0.5(1 ± 77 ) appears due to 
the application of 77oe in the definition of Figs. ([2]) 
and © demonstrate that the soliton width A is of order 
of lO^rjje = 5 X 10“^C772. Increasing the soliton velocity 
Uo we decrease the soliton width A. Since the soliton 
width A should be larger then the average interparticle 
distance A ^ we obtain that Uq 10~‘^VFe and 

77is ^ 10^^ cm“^. 

Amplitude of the dark soliton in the subsystem of spin- 
down electrons increases with the decrease of the spin 
polarization Amplitude of the bright soliton in the 
subsystem of spin-up electrons has nonmonotonic depen¬ 
dence on the spin polarization ([3]). It has the minimum 
at the intermediate spin polarization 770 = 0.54. The 
amplitude increases at change of the spin polarization 
77 from 770 towards smaller or larger values. At change 
of the spin polarization from 770 = 0.54 to 77 = 0.1 and 
77 = 0.9 the amplitude increases from 4 to 6.7 and 8.2 
correspondingly. 

We present analysis of the soliton characteristics in the 
self-consistent field approximation. In the next section 
we include the contribution of the Coulomb exchange in¬ 
teraction in the spin-electron acoustic soliton propaga¬ 
tion. 


V. CONTRIBUTION OF EXCHANGE 
INTERACTION IN SPIN ELECTRON ACOUSTIC 
SOLITON 


The spin dependence of the Coulomb exchange inter¬ 
action force field in the electron gas is calculated in Ref. 
[ 3 ^ . Significant role of the Coulomb exchange interaction 
at rather large spin polarization is demonstrated there. 

To underline the fact that, in this section, we consider 
the Coulomb interaction beyond the self-consistent field 
approximation we present the explicit form of the per¬ 
turbation velocity obtained in Sec. Ill (see formula (j35j)) 
= no-^Uf + noiU^/noe, where 77oe = ^ot + func¬ 
tions and are defined by formulae (|32]) and (|33]) . 
correspondingly. 

Coefficient A of the KdV equation (EZl) at the account 
of the Coulomb exchange interaction appears as 

^ ^ ^ - ul)^ }■ 

The exchange interaction between spin-down electrons 
being in states occupied by single electron modifies Uy 
The account of the Coulomb exchange interaction does 
not change sign of A, so we have A > 0. 
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The contribution of the exchange interaction in coeffi¬ 
cient 
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(y2 _ ^2)3 6 me J 


(43) 


reveals in several ways. We find modifications of (see 
the second term in formula (|33|) i and V‘^ (compare for¬ 
mulae §B)- Moreover, we find an extra term (the last 
term in formula (|43|) ) in coefficient B. 

Coefficient B is positive in the self-consistent field ap¬ 
proximation Bscf > 0 . The dependence of the Bscf on 
the electron concentration is located in Bq. If we include 
the exchange interaction the behavior of the coefficient 
B becomes rather complicate. We find a dependence of 
B/Bq on the equilibrium electron concentration noe- 

This dependence reveals in properties of the soliton 
profile described below. We can track this dependence at 
an intermediate spin polarization r]' = 0.5. We see that 
coefficient B is positive B > 0 in the regime of rather 
large concentrations noe ^ 10 ^^ cm“^, which is simi¬ 
lar to the self-consistent field approximation. At ho ~ 
5.4 X 10^^ cm“^ we find fast increase of B {B ^ +oc). 
After point ho the coefficient B increases from minus in¬ 
finity up to the zero value 5 = 0 at Uq = 0.94 x 10^^ 
cm“^. At smaller concentrations the coefficient B is pos¬ 
itive. We see that the coefficient 5 as a function of the 
electron concentration demonstrates the hyperbolic de¬ 
pendence. Since B is in the denominator of the soli¬ 
ton amplitude, the soliton amplitude tents to zero at ho- 
Hence, the soliton does not exist near this concentra¬ 
tion. At noe = tiq coefficient B vanishes. Consequently, 
the soliton amplitude approach infinity. Therefore, the 
perturbation method, we apply to find the spin-electron 
acoustic soliton, cannot give any information about soli¬ 
ton behavior near this point. 

For the formation of the soliton we need to have waves 
with the stable linear spectrum. Hence, we need to have 
a positive square of the perturbation velocity given by 
formula (|35]) . which is a combination of and U^. In 
the self-consistent field approximation considered in the 
previous section the spectrum is stable for all values of 
parameters. If we include the exchange interaction sit¬ 
uation changes. Parameter contains a negative term 
caused by the exchange interaction (see formula (|33]) b 
Therefore, V‘^ can become negative. To find areas of 
positive we present Fig. [5l In the left-hand column 
we depict and as the functions of the spin 

polarization 77 at the different equilibrium electron con¬ 
centrations noe- The lower row of pictures is obtained 
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FIG. 5: (Color online) The left column on the figure shows 
the dependencies of Wq = /'^‘pe (thick black solid line), 

Wi = 3F|(?7)/7;|^e (thin red solid line), W 2 — ‘^Ul(inj) 

(blue dashed line) at different parameters of the electron gas 
(in order from the upper picture to the lower picture we have 
no = 10^^ cm“^ with the exchange interaction, no = 10^^ 
cm“^ with the exchange interaction, no = 10^^ cm“^ with the 
exchange interaction, any no with no account of the exchange 
interaction). In this figure Wa stands for {Wq , VFf, Wl}• The 
right column on the figure shows differences of the velocity 
squares presented in denominators in formulae ( 1301 ) and m 
VFo - VFf (red solid line) and VFo^ - VF| (blue dashed line). 
In the figure we apply notation Ws for Wi and . Regimes 
for differences Wq — Wg in the right column correspond to 
the regimes for in the left column. 


in the self-consistent field approximation to provide the 
comparison with the results of the previous section. 

The results of the previous section show that sign of 
quantities V‘^ — and V‘^ — define profiles of the 
spin-electron acoustic solitons of the concentration of the 
spin-up and spin-down electrons. Therefore, we present 
these quantities in the right-hand column in Fig. O 

Fig. [5] shows that at noe = 10^^ cm“^ results includ¬ 
ing the exchange interaction have small, but noticeable, 
difference with the results of the self-consistent field ap¬ 
proximation. Hence, in this regime, we find the soliton 
profiles (see the lower row in Fig. E]) similar to the pro¬ 
files found in the previous section (see Figs. [2]and[3|). 
However, we see the contribution of spin-down electrons 
presented by blue dashed lines in Fig. [5] on 14 percent. 
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FIG. 6: (Color online) The figure shows the soliton profiles 
for spin-down (the left column) and spin-up (the right col¬ 
umn) electron concentrations for different equilibrium con¬ 
centrations of electrons presented in the figure. We present 
the soliton profile in terms of = (nis/nos)/{V^Uo/vFe), 
where nis appears at the substitution of solution (|38l) in for¬ 
mulae m and (|30)). 


approximately, due to the account of the exchange inter¬ 
action. It reveals in the amplitude of the spin-down soli¬ 
ton. Linear spectrum is stable in this regime > 0^ and 
the parameters V‘^ — and V‘^ — have positive and 
negative signs correspondingly {U^ < < U^). In this 

regime the coefficient B is positive B > 0. Consequently, 
we find the dark soliton in the spin-down electron con¬ 
centration and the bright soliton in the spin-up electron 
concentration. 

Decreasing the equilibrium concentration of electrons 
down to noe = 10^^ cm“^ we find that the spectrum is 
stable {V‘^ > 0) for all 77 G [0,1]. However, relative val¬ 
ues of are changed. Velocity Ui becomes the 

smallest of them. So, we have < U^. Conse¬ 

quently the signs of parameters V‘^ — and V‘^ — are 
changed either. In this regime the coefficient B is neg¬ 
ative B < 0. Therefore, we find the dark soliton in the 
spin-down electron concentration and the bright soliton 
in the spin-up electron concentration. 

At noe = 10^^ cm“^ the linear spectrum becomes un¬ 
stable in the wide range of the spin polarization. We find 
the stability interval at 77 G (0.01,0.06). In this regime 
the coefficient B is positive B > 0. Consequently, in the 
area of stability, we find the bright soliton in the spin- 
down electron concentration and the dark soliton in the 
spin-up electron concentration (see the upper row in Fig. 


El). 


VI. CONCLUSIONS 

Existence of the spin-electron acoustic soliton has been 
discovered. Its existence closely related to the spin elec¬ 
tron acoustic waves recently obtained in Ref. The 

SEAW has linear spectrum, so it resembles similarity to 
the ion-acoustic wave, but the SEAW has larger frequen¬ 
cies. The balance between the dispersion and nonlin¬ 
earity in the SEAWs of small amplitude allows to form 
a soliton solution obtained in this paper and called the 
spin-electron acoustic soliton. 

To study the non-linear spin-electron acoustic waves 
we have developed a generalization of the separate spin 
evolution quantum hydrodynamics. This generalization 
includes the Coulomb exchange interaction. The ex¬ 
change interaction appears from the interaction of the 
spin-down electrons being in the states occupied by one 
electron only. This mechanism was considered in the 
single fluid model of electrons. In this paper we have 
adopted it for the SSE-QHD. 

We have considered the spin-electron acoustic soliton 
in two regimes. Eirst of all we have considered it in, 
rather simple, regime of the self-consistent field approx¬ 
imation. We have found that the spin-electron acoustic 
soliton shows itself as the dark soliton of the spin-down 
electron concentration and the bright soliton of the spin- 
up electron concentration. The soliton shows similar be¬ 
havior for all equilibrium concentrations of electrons. 

The second regime of the spin-electron acoustic wave 
study includes the Coulomb exchange interaction. 

The exchange interaction significantly change proper¬ 
ties of the spin-electron acoustic soliton. Strong depen¬ 
dence of the soliton properties reveals in this regime. 
At equilibrium concentration noe = 10 ^^ cm“^ we have 
found, in opposite to the self-consistent field approxima¬ 
tion, the bright soliton of the spin-down electron con¬ 
centration and the dark soliton of the spin-up electron 
concentration, existing in the narrow interval at rather 
small spin polarizations 77 G (0.01,0.06). 

The increase of the equilibrium electron concentration 
increases the area of the soliton existence. At 77oe = 10^^ 
cm“^ we have found the existance of the spin-electron 
acoustic soliton at 77 G (0.01,0.99). At 77oe > 10^^ cm“^ 
we have obtained the dark soliton of the spin-down elec¬ 
tron concentration and the bright soliton of the spin-up 
electron concentration. It is in the agreement with the 
self-consistent field approximation. However, the param¬ 
eters of the soliton at the account of the exchange inter¬ 
action at noe ^ [ 10 ^^, 10 ^^] cm“^ differ from the results 
of the self-consistent field approximation. 
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